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VOLTAGE-CONTROLLED  FERROELECTRIC  LENS  PHASED  ARRAYS 


1.  INTRODUCTION 

Phased  array  antennas  can  steer  transmitted  and  received  signals  without  mechanically  rotating  the 
antenna.  Each  radiating  element  of  a  phased  array  is  normally  connected  to  a  phase  shifter  and  a  driver, 
which  determines  the  phase  of  the  signal  at  each  element  to  form  a  beam  at  the  desired  angle.  The  most 
commonly  used  phase  shifters  are  ferrite  and  diode  phase  shifters.  Phase  shifters  using  ferroelectric 
materials  have  been  proposed  [1,2].  The  cost  of  a  phased  array  mainly  depends  on  the  cost  of  phase 
shifters  and  drivers.  A  typical  array  may  have  several  thousand  elements  and  that  many  phase  shifters  and 
drivers;  hence,  it  is  very  expensive.  Therefore,  reducing  the  cost  and  complexity  of  the  phase  shifters, 
drivers,  and  controls  is  an  important  consideration  in  the  design  of  phased  arrays.  The  phased  array 
reported  here  uniquely  incorporates  bulk  phase  shifting;  the  array  does  not  contain  individual  phase 
shifters  but  uses  ferroelectric  material.  This  will  reduce  the  number  of  phase  shifters  from  (n  x  m)  to 
(n  +  m),  where  n  is  the  number  of  columns  and  m  is  the  number  of  rows  in  a  phased  array.  The  number  of 
phase  shifter  drivers  and  phase  shifter  controls  is  also  significantly  reduced  by  using  row-column  beam 
steering.  Bulk  phase  shifting  using  diodes  has  been  proposed  [3]  and  reasonably  developed  [4]  in  the 
Radant  lens.  The  lens  described  here  uses  voltage-controlled  ferroelectric,  which  introduces  analog  phase 
shift  rather  than  digital  phase  shift  as  in  the  Radant  lens.  The  ferroelectric  lens  has  the  further  advantages 
of  smaller  lens  thickness,  higher-power  handling  capability,  and  simpler  beam  steering  controls,  and  it 
uses  less  power  to  control  the  phase  shift  compared  to  the  Radant  lens.  Thus,  it  leads  to  low-cost  phased 
arrays.  However,  it  should  be  noted  that  the  use  of  row-column  steering  may  limit  the  level  of  side  lobes 
that  can  be  achieved.  For  the  present  study,  it  is  assumed  that  ultralow  side  lobes  are  not  a  requirement 
and,  therefore,  row-column  phase  control  can  be  used  to  significantly  reduce  the  phased  array  cost. 

2.  DESCRIPTION  OF  FERROELECTRIC  LENS  AND  ITS  OPERATION 

The  main  feature  of  the  antennas  that  use  ferroelectric  material  is  the  change  of  permittivity  with  an 
applied  DC  control  voltage.  A  lens-type  antenna  is  discussed  in  this  report.  Figure  1  shows  a  dielectric 
lens  made  of  dielectric  slabs  sandwiched  between  conducting  plates.  Dielectric  slabs  are  made  of 
ferroelectric  material  whose  dielectric  constant  can  be  changed  by  applying  and  varying  the  DC  electric 

field  (DC  voltage  sources  Vj,  V2 . V„  are  used  for  this  purpose,  as  shown  in  Fig.  1).  If  a  plane  wave  is 

incident  on  one  side  of  the  lens  with  RF  electric  field  E  normal  to  the  conducting  plates,  the  beam  coming 
out  on  the  other  side  of  the  lens  can  be  scanned  in  the  E-plane  if  a  linear  phase  gradient  is  introduced 

along  the  E-plane  direction  by  adjusting  the  voltages  V|,  V2 . V„.  The  corresponding  dielectric  constants 

are  shown  as  ej,  £2, ....£„  in  Fig.  1. 

Figure  1  illustrates  the  principle  of  operation  of  a  ferroelectric  lens.  However,  it  has  some  limit¬ 
ations.  To  reduce  the  lens  complexity  and  weight,  the  number  of  DC  control  voltages  and  the  number  of 
conducting  plates  must  be  minimized.  This  can  be  achieved  by  limiting  the  space  between  the 
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conducting  plates  to  be  less  than  Afl/(1  +  sin  (t)^),  similar  to  inter-element  spacing  in  any  phased  array,  to 
avoid  grating  lobes,  where  Xq  is  the  free-space  wavelength,  and  (|)s  is  the  maximum  scan  angle.  Normally, 
the  space  between  the  conducting  plates  would  be  -Xq/I.  However,  in  a  ferroelectric  lens,  this  space  is 
filled  with  a  high-dielectric  constant  ferroelectric  material,  and  the  space  between  the  plates  is  much 
larger  than  X/2,  where  X  is  the  wavelength  in  the  ferroelectric.  This  means  that  higher-order  modes  may 
propagate.  Therefore,  to  eliminate  the  problem  with  higher-order  modes,  the  thickness  of  the  ferroelectric 
material  should  be  kept  to  less  than  X/2.  In  addition,  there  should  be  some  type  of  impedance  matching 
arrangement  to  match  the  lens  surface  to  free  space. 

Figure  2  shows  a  more  practical  lens  configuration.  This  configuration  has  several  advantages  over 
the  basic  configuration  of  Fig.  1.  A  stepped  configuration  is  used  for  impedance  matching  to  reduce  the 
amount  of  dielectric  material,  to  eliminate  higher-order  mode  propagation,  and  to  reduce  the  bias  voltage 
necessary  to  create  a  certain  amount  of  DC  electric  field  intensity  within  the  ferroelectric  material. 
Quarter-wave  transformers  are  used  for  impedance  matching  the  empty  waveguide  region  to  the 
ferroelectric  loaded  waveguide  region.  More  than  one  transformer  may  be  needed,  depending  on  the  value 
of  the  dielectric  constant  of  the  ferroelectric  material  and  the  degree  of  matching  needed.  Figure  2  also 
shows  a  further  refinement.  This  refers  to  the  way  the  DC  bias  voltage  is  applied  to  the  ferroelectric 
material.  An  additional  center  conducting  plate  is  used  to  bifurcate  the  ferroelectric  material.  The  DC 
voltage  is  applied  between  this  center  conducting  plate  (recessed  with  respect  to  parallel  plates)  and  the 
parallel  plates.  Since  the  dielectric  material  is  bifurcated,  only  half  the  DC  voltage  needs  to  be  applied  to 
produce  the  same  DC  electric  field  intensity.  It  is  this  DC  electric  field  intensity  that  controls  the 
dielectric  constant  of  the  ferroelectric.  In  addition,  the  parallel  plates  are  at  ground  potential.  This  makes 
the  handling  of  the  lens  safer. 
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Parallel  Conducting  Plates 


(a)  Ferroelectric  lens 


(b)  Top  view  of  one  column 


Fig.  2  —  A  more  practical  configuration  of  the  ferroelectric  lens 

The  lens  thickness  needed  to  use  the  lens  as  a  scanning  antenna  can  be  determined  as  follows:  let  t  be 
the  lens  thickness  as  shown  in  Fig.  1,  and  let  the  dielectric  constant  of  any  one  dielectric  slab  be 
when  no  bias  voltage  is  applied.  The  dielectric  constant  of  a  ferroelectric  decreases  as  the  bias  voltage  is 
increased.  Let  ^nin  dielectric  constant  when  maximum  DC  bias  is  applied.  The  phase  introduced 

by  the  slab  when  no  bias  voltage  is  applied  is  given  by 
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Similarly,  the  phase  introduced  by  the  same  slab,  when  maximum  bias  voltage  is  applied,  is  given  by 

<l’min  =T^^V^r,imn  •  (2) 

Aq 

Note  the  dielectric  constant  decreases  with  applied  bias  voltage.  Then  the  maximum  phase  change  when 
the  maximum  bias  voltage  is  applied  is 

^max  ~  V^r,niin  )* 

Aq 

For  scanning  applications,  this  maximum  phase  change  should  be  equal  to  2tc  (as  in  the  case  of  a  phase 
shifter  used  in  a  scanning  array);  hence,  Eq.  (3)  becomes 

2^1  =  ^(^^r,max  "■  ^j^r, min  j  • 


So,  the  thickness  of  the  ferroelectric  material  needed,  in  the  direction  of  propagation,  to  obtain  360°  phase 
shift  is 

t  =  ,  ^0,  -  .  (5) 


^j^r, max 


Tunability  can  be  defined  for  a  ferroelectric  as  the  fractional  change  in  the  dielectric  constant  with  applied 
DC  bias  voltage  or 


Tunability  = 


^A^max  ^r,min 


Now,  Eq.  (5)  can  be  written  in  terms  of  tunability  as 


^0  V^r,max[l-Vl-t“”abilty] 


So,  the  thickness  of  the  ferroelectric  material  needed  is  a  function  of  the  dielectric  constant  and  the 
tunability  of  the  ferroelectric  and  the  wavelength.  The  relationship  shown  in  Eq.  (7)  is  illustrated  in 
Fig.  3  for  typical  of  90  and  120.  Additional  discussion  on  this  relationship  is  included  later. 

Also,  it  can  be  shown  [5]  that  the  attenuation  constant  for  a  uniform  plane  wave  propagating  in  a  low- 
loss  dielectric  is 


a  (dB)  =  27.3  tan  6 


where  tan  6  is  the  loss  tangent.  Using  Eq.  (7),  we  see  that  to  obtain  360”  phase  shift,  the  dielectric  loss 
through  the  ferroelectric  is 


a(dB)  = 


27.3  tang 
1  -  Jl  —  tunability 


(9a) 
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Fig.  3  —  Amount  of  ferroelectric  needed  for  360°  phase  shift 


It  may  be  noted  that  the  lens  loss  is  independent  of  the  ferroelectric  permittivity  and  depends  only  on  its 
loss  tangent  and  tunability.  Figure  4  plots  this  equation.  It  can  also  be  approximated  by  using  the  binomial 
series  as 


a(dB)  = 


55  tan  5 
tunability 


(9b) 


Fig.  4  —  Ferroelectric  loss  for  360°  phase  shift 
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In  general,  the  ferroelectrics  with  higher  dielectric  constant  offer  higher  tunability,  which  is  desired  to 
reduce  the  lens  thickness.  However,  matching  the  lens  to  free  space  is  easier  for  smaller  8,..  Therefore,  a 
compromise  is  needed  between  reducing  the  lens  thickness  (to  reduce  overall  lens  size)  and  achieving 
reasonable  impedance  match  (to  reduce  reflections  from  the  lens  surface).  For  a  typical  value  of  ~  100, 
it  is  possible  to  obtain  a  tunability  of  20  to  30%,  which  results  in  a  reasonable  lens  thickness  of  ~  (e.g., 

3  cm  at  10  GHz),  as  shown  in  Fig.  3.  For  this  reason,  ferroelectric  materials  with  8^  ~  100  were  chosen  for 
the  experimental  work,  which  will  be  discussed  later.  Figure  4  shows  that  a  tan  8  =  0.005  is  needed  to 
limit  the  lens  loss  to  less  than  1  dB.  The  existing  ferroelectric  materials  are  a  bit  more  lossy  (tan  5  =  0.008 
at  10  GHz)  and  further  research  is  needed  to  reduce  losses.  Additional  discussion  regarding  this  will  be 
included  when  the  experimental  results  are  discussed  in  a  later  part  of  the  report. 

The  above  theory  can  be  applied  to  any  dielectric  whose  dielectric  constant  can  be  changed  by  an 
external  stimulus.  In  addition  to  ferroelectric  materials  (which  are  the  focus  of  this  report),  liquid  crystals 
are  also  voltage  variable  dielectric  materials.  The  dielectric  constant  of  a  liquid  crystal  can  be  changed  by 
applying  a  bias  voltage,  and  thus  liquid  crystals  can  also  be  used  as  phase  shifters  [6].  Similar  to  the 
tunability  of  a  ferroelectric,  a  liquid  crystal  has  a  birefringence  defined  as 

~  ■\j^r,max  ~  •  (1^^) 

Unlike  the  ferroelectric,  however,  the  dielectric  constant  of  the  liquid  crystal  increases  with  bias  voltage. 
So,  8^  is  the  dielectric  constant  of  the  liquid  crystal  when  maximum  DC  bias  is  applied,  and  8^ 
is  the  dielectric  constant  when  no  bias  voltage  is  applied.  Therefore,  the  thickness  of  the  liquid  crystal 
needed  (in  the  direction  of  propagation)  to  obtain  360°  phase  shift  is 


and  the  dielectric  loss  through  the  liquid  crystal  to  obtain  360°  phase  shift  is 


a(dB)  = 


27.3  tan  S-sJ^ 
An 


(10c) 


It  is  not  the  purpose  of  this  report  to  present  the  application  of  liquid  crystals  to  phased  array 
antennas.  Nevertheless,  this  report  will  provide  considerable  insight  into  the  design  of  an  electronically 
scanning  lens  antenna  based  on  liquid  crystals  because  the  basic  principles  and  design  equations  are 
identical  to  those  of  a  ferroelectric  lens.  We  will  end  our  discussion  of  liquid  crystals  by  stating  that  the 
main  disadvantage  of  liquid  crystals  at  microwave  frequencies  is  the  slow  response  time  of  the  material 
when  the  bias  is  removed. 

The  ferroelectric  lens  offers  electronic  scanning  in  one  plane.  The  lens  proposed  here  can  be  fed  by  a 
nonscanning  planar  array,  like  a  slotted  waveguide  array.  A  combination  of  slotted  waveguide  array  with 
phase  shifters  and  the  lens  proposed  here  can  be  used  as  a  phased  array  that  can  scan  in  two  planes.  A 
space  feed  can  be  used  with  the  combination  of  two  lenses  proposed  here  (with  a  polarization  rotator  in 
between)  to  scan  the  beam  in  two  planes.  The  details  of  these  alternatives  are  discussed  in  the  next 
section. 
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3.  PHASED  ARRAY  CONFIGURATIONS  USING  FERROELECTRIC  LENS 

FOR  TWO-  DIMENSIONAL  SCANNING 

A  single  ferroelectric  lens  provides  one-dimensional  (1-D)  electronic  scanning.  Two-dimensional 
(2-D)  scanning  can  be  achieved  by  cascading  two  ferroelectric  lenses  or  using  one  ferroelectric  lens  in  a 
hybrid  configuration  with  a  phased  array  that  can  scan  in  one  plane. 

Figure  5  shows  the  cascading  of  two  spatially  orthogonal  ferroelectric  lenses.  The  first  lens  provides 
an  elevation  scan  of  a  vertically  polarized  wave.  A  90°  passive  polarization  rotator  then  rotates  the 
electric  field  to  become  horizontally  polarized.  The  second  lens  then  provides  the  azimuth  scanning  of  the 
horizontally  polarized  wave.  In  Fig.  5,  a  planar  array  is  shown  as  the  illuminator  (or  feed)  for  the  dual 
lenses.  A  space  feed  can  also  be  used,  as  shown  in  Fig.  6,  with  the  dual  lens  configuration.  In  this 
configuration,  in  addition  to  scanning  the  beam,  row-column  phase  controls  can  also  be  used  to  correct 
the  spherical  phase  errors  due  to  the  point  space  feed;  however,  this  phase  correction  is  not  exact,  but  it  is 
satisfactory  for  many  applications.  Spherical  phase  error  corrections  for  an  offset  feed  configuration  is  not 
possible  with  row-column  phase  controls.  So,  offset  feeding  is  not  recommended  for  dual  lens 
configuration. 


Polarization  Rotator 


Fig.  5  —  Dual  lens  configuration  with  a  planar  array  feed 


Another  method  of  achieving  2-D  scanning  uses  a  hybrid  technique  in  which  a  conventional  planar 
array  with  phase  shifters  provides  scanning  in  one  plane  and  the  ferroelectric  lens  provides  the  scanning 
in  the  other  plane.  For  example,  consider  a  stacked  slotted  waveguide  array  with  the  slots  in  the  narrow 
wall  of  the  waveguide.  Each  waveguide  is  fed  through  a  phase  shifter  at  its  input.  Scanning  in  one  plane 
is  achieved  with  these  phase  shifters.  Scanning  in  the  orthogonal  plane  is  obtained  by  placing  a 
ferroelectric  lens  in  front  of  the  waveguide  array,  as  shown  in  Fig.  7.  For  the  reasons  mentioned  before, 
the  cost  of  the  phased  arrays  that  use  ferroelectric  lens  will  be  much  less  than  the  phased  arrays  that  use  a 
phase  shifter  or  a  transmit/receive  (T/R)  module  behind  each  radiating  element. 
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Lens  #2 


Lens  #1 
Feed  Horn 


Fig.  6  —  Dual  lens  configuration  with  a  space  feed 
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Fig.  7  —  Hybrid  phased  array  configuration  for  2-D  scanning 
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4.  FERROELECTRIC  MATERIALS 

Ferroelectric  materials  are  nonlinear  dielectrics.  They  possess  spontaneous  polarization  similar  to 
ferromagnetic  materials.  However,  in  our  application,  we  use  the  material  in  the  nonferroelectric 
(paraelectric)  state  above  the  Curie  temperature,  where  it  is  isotropic  and  provides  reciprocal  phase  shift. 
The  high  and  tan  6  of  pure  ferroelectric  materials  have  limited  their  applications  in  phased  arrays. 
Recent  improvements  in  ferroelectric  composite  materials  at  the  Army  Research  Laboratory  (ARL), 
Aberdeen  Proving  Ground,  Maryland,  prompted  renewed  interest  in  developing  phase  shifters  at 
microwave  frequencies.  The  materials  are  bulk  oxide-ceramic  composites  of  barium  strontium  titanate 
oxide  (BSTO).  ARL  has  also  developed  a  simple,  low-cost  processing  technique  for  their  manufacturing 
[7-9].  It  involves  mixing  barium  titanate  (BaTi03)  and  strontium  titanate  (SrTi03)  powders.  The  resulting 
BSTO  powder  is  calcined  and  mixed  with  an  oxide.  After  pressing  and  sintering,  the  composite  is  ready 
to  be  machined.  Table  1  lists  the  properties  of  several  ARL  ferroelectric  composites  at  room  temperature. 
The  first  column  indicates  the  Ba/Sr  ratio  of  the  composition.  The  second  column  shows  the  weight 
percent  oxide  in  the  composition.  The  third  column  is  the  Curie  temperature.  The  fourth  column  is  the 
ferroelectric  tunability  at  2  V/jim.  The  lower  the  Curie  temperature,  the  more  paraelectric  and  less 
tunable  the  material  becomes  at  room  temperature.  The  last  two  columns  list  and  tan  6  measured  at  1 
kHz  and  10  GHz.  Low-frequency  results  were  calculated  from  capacitance  measurements.  Measurements 
at  10  GHz  were  performed  by  Dr.  R.  Geyer  at  the  National  Institute  of  Standards  and  Technology  [10] 
and  confirmed  with  our  measurements.  We  see  that,  in  general,  tunability,  8^.,  and  tan  5  decrease  with 
decreasing  Ba  content;  they  also  decrease  with  increasing  oxide  content.  After  oxide  is  added,  changes 
much  less  with  firequency.  As  expected  for  ceramics,  tan  8  increases  with  frequency;  however,  it  is  far 
lower  than  that  of  any  other  tunable  ferroelectrics.  Measurements  at  higher  bias  voltages  indicate  that 
tunability  increases  linearly  with  an  increase  in  bias  voltage. 


Table  1  —  Properties  of  ARL  Bulk  Ferroelectrics 


Material 

Oxide 
(wt  %) 

Curie 

Temp. 

Tunability 

(%) 

8, 

1  kHz/ 10  GHz 

tan  5x  10~3 

1  kHz  /  10  GHz 

®^.5oS’‘o.5o"^^3 

0 

-25°C 

25.5 

1908/1099 

20.0/18.9 

20 

-55°C 

6.44 

592/616 

0.73/8.7 

1  B^.5oSro.5oTi3 

30 

-70°C 

5.82 

414/463 

0.7/8.4 

60 

-95°C 

3.66 

78/84 

0.76/6.5 

®3o.55Sro.45Ti3 

30 

-45°C 

9.54 

478/527 

0.75/12.1 

®^.55S%45Ti3 

60 

-50°C 

6.46 

95/100 

0.34/7.9 

60 

-55°C 

9.99 

117/118 

1.48/12.9 

5.  THEORETICAL  ANALYSIS  AND  EXPERIMENTAL  RESULTS 

For  the  theoretical  analysis  of  the  ferroelectric  lens,  an  individual  section  between  two  conducting 
parallel  plates  of  the  lens  can  be  considered  as  one  column  of  a  phased  array.  Initially,  the  column  height 
is  considered  infinite  in  extent  and  hence  it  becomes  a  2-D  parallel  plate  waveguide  with  electric  field 
normal  to  the  plates  as  shown  in  Fig.  2(b).  The  lens  column  fed  by  a  parallel  plate  waveguide  was 
evaluated  assuming  transverse  electromagnetic  (TEM)  mode  propagation.  To  eliminate  higher-order 
mode  propagation,  the  ferroelectric  slab  thickness,  between  one  of  the  parallel  plates  and  the  central 
conducting  plate,  is  limited  to  less  than  XU,  as  discussed  earlier  and  shown  in  Fig.  2(b).  Theoretical 
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analysis  was  also  performed  using  Touchstone,  a  commonly  available  software,  for  a  rectangular 
waveguide  with  a  cross  section  shown  in  Fig.  2(b).  The  analysis  was  done  assuming  different  column 
heights  (waveguide  widths).  For  widths  larger  than  4Xo,  the  rectangular  waveguide  results  matched  well 
with  parallel  plate  results.  For  practical  reasons,  most  of  the  initial  measurements  were  performed  in  a 
reduced  height  X-band  rectangular  (WR90)  waveguide.  Theoretical  and  experimental  results  for  the 
waveguide  configuration,  along  with  a  brief  description  of  the  ferroelectric  lens  and  its  principle  of 
operation,  have  been  reported  earlier  [11-14].  This  report  includes  the  detailed  theoretical  analysis  of  the 
ferroelectric  lens,  along  with  the  extensive  experimental  results.  Theoretical  analysis  also  indicated  that 
the  results  are  better  for  a  parallel  plate  waveguide  than  for  a  rectangular  waveguide. 

In  this  report,  for  theoretical  analysis  and  numerical  results,  the  geometry  was  restricted  to  be  2-D 
(planar)  since  each  column  of  the  ferroelectric  lens  is  predominantly  2-D,  as  discussed  earlier.  Because 
the  lens  is  to  be  fed  with  a  linearly  polarized  field  with  the  electric  field  normal  to  the  conducting  plates, 
the  lens  was  analyzed  only  for  that  case.  A  single  column  of  the  lens,  fed  by  a  parallel  plate  waveguide 
was  evaluated  as  a  parallel  plate  waveguide.  Therefore,  the  results  presented  here  are  for  a  single  column 
of  the  lens  and  do  not  include  the  effects  of  mutual  coupling  between  columns. 


Analytical  Solution 


We  have  considered  two  impedance  matching  networks  to  match  the  lens  surface  to  free  space  over  a 
broad  range  of  frequencies.  Both  networks  use  multisection  impedance  matching  transformers.  They  are 
shown  in  Fig.  8,  which  are  the  top  views  of  one  column.  The  matching  network  in  Fig.  8(a)  uses  three 
transformers.  The  matching  network  in  Fig.  8(b)  uses  four  transformers.  Each  transformer  is  a  section  of  a 
parallel  plate  waveguide.  In  Fig.  8(a): 


where 


In  Figure  8(b): 


where 


Li  =L^_i,bi  =  b^_i,Zi  -88_,-,e4  -£^,^2  -  b-i  -b^  -  b, 

1  =  1,2,  3. 

1  =  1,2,  3,4. 


Li  are  the  lengths  of  the  waveguide  sections;  bi  are  the  heights  (separations  between  the  parallel 
plates)  of  the  waveguide  sections;  are  the  complex  permittivities  of  the  dielectrics  completely  filling  the 
waveguide  sections;  and  is  the  dielectric  constant  of  the  ferroelectric  material.  Since  the  input  and 
output  sections  of  the  lens  are  identical  (empty  parallel  plate  waveguide),  the  same  matching  network  can 
be  used  at  both  ends.  The  matching  network  provides  impedance  matching  between  the  empty  semi¬ 
infinite  parallel  plate  waveguide  (A  and  C  in  Fig.  8)  and  the  waveguide  section  completely  filled  with  the 
ferroelectric  (B  in  Fig.  8).  The  heights  of  waveguides  A,  B,  and  C  are  chosen  to  assure  that  only  the 
dominant  TEM-mode  will  propagate  in  the  waveguides  and  that  grating  lobes  are  not  formed  when  an 
antenna  is  designed  by  arraying  several  columns  together.  Because  waveguide  B  is  filled  with  a  high- 
dielectric-constant  material,  its  height  is  chosen  to  be  much  smaller  than  the  height  of  waveguides  A  and 
C  to  eliminate  higher-order  mode  problems.  The  heights  of  waveguides  A  and  C  are  the  same  and  are  less 
than  Xq/2  at  the  highest  frequency  of  operation  in  order  to  avoid  grating  lobes. 
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I  II  III  IV  V  VI  vii  vm 


ABC 

(a) 


I  II  III  IV  V  VI  VII  VIII  IX  X 


J  i . . . 1  i . - . . . 

ABC 

(b) 


Li  are  lengths  of  waveguide  sections 
hi  are  the  waveguide  heights 


Ax 


z 


€i  are  the  dielectric  constants 

Fig.  8  —  Matching  networks:  (a)  three-section  matching  network, 
and  (b)  four-section  matching  network 


The  matching  network  of  Fig.  8(a)  uses  a  design  with  three  transformers,  where  the  first  transformer 
is  an  empty  parallel  plate  waveguide  whose  height  is  between  that  of  waveguide  A  (or  C)  and  B.  The 
other  two  transformers  are  parallel  plate  waveguides  completely  filled  with  dielectric  materials  whose 
relative  permittivities  (e^)  are  between  that  of  air  (or  unity)  and  that  of  the  ferroelectric.  Dielectric 
constant  of  the  first  transformer  is  less  than  that  of  the  second  transformer.  The  heights  of  these  two 
transformers  are  the  same  as  the  height  of  waveguide  B.  The  matching  network  of  Fig.  8(b)  uses  four 
transformers,  where  the  first  two  transformers  are  empty  parallel  plate  waveguides  whose  heights  are 
between  that  of  waveguide  A  (or  C)  and  B.  The  height  of  the  first  transformer  is  larger  than  the  height  of 
the  second  transformer.  The  other  two  transformers  are  parallel  plate  waveguides  completely  filled  with 
dielectric  materials  whose  relative  permittivities  are  between  that  of  air  (or  unity)  and  that  of  the 
ferroelectric.  The  heights  of  these  two  transformers  are  the  same  as  the  height  of  waveguide  B.  Dielectric 
constant  of  the  first  transformer  is  less  than  that  of  the  second  transformer. 


Since  a  single  ferroelectric  lens  provides  electronic  scanning  in  one  plane,  another  lens  or  antenna  is 
needed  to  scan  in  the  second  plane.  If  the  first  antenna  (the  one  closer  to  the  receiver/transmitter)  is  a  lens, 
the  electromagnetic  wave  is  incident  normally  on  this  lens  from  the  feed.  This  lens  will  create  a  phase 
gradient  in  one  plane,  and  its  radiation  will  be  incident  obliquely  onto  the  second  lens.  For  the  theoretical 
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analysis  of  one  column,  we  assume  that  the  dominant  TEM-mode  is  obliquely  incident  in  the  ^T-plane. 
The  electric  field  is  orthogonal  to  the  waveguide  plates  (y-direction).  The  magnetic  field  is  in  the  plane 
parallel  to  the  waveguide  walls.  The  oblique  incident  angle  is  measured  in  the  x-z  plane.  This  is  the  scan 
angle  of  the  first  lens.  A  time  dependency  of  expO'coO  is  assumed,  and  the  waveguide  walls  are  assumed 
to  be  perfect  conductors. 

As  a  plane  TEM  wave  is  obliquely  incident  on  the  matching  network  from  waveguide  A,  some  of  the 
energy  associated  with  the  incident  wave  will  be  reflected  back  to  waveguide  A.  We  will  derive  an 
analytical  solution  for  the  reflection  coefficient  and  use  this  solution  to  design  the  optimum  matching 
network,  which  minimizes  the  magnitude  of  the  reflected  field  for  a  range  of  incidence  angles  over  the 
frequency  band  of  interest  and  for  the  required  change  of  the  ferroelectric  dielectric  constant  to  obtain 
360°  differential  phase  shift.  The  case  of  normal  incidence  represents  a  special  case  of  oblique  incidence 
where  the  incidence  angle  is  zero. 

The  optimization  algorithms  were  devised  and  computer  programs  have  been  written  based  on  the 
analytical  solutions.  As  optimization  parameters,  computer  programs  use  lengths  of  the  transformers, 
height(s)  of  the  empty  waveguide  transformer(s),  and  the  dielectric  constants  of  the  materials  in  the  last 
two  transformers. 


First,  consider  the  matching  network  which  uses  four  transformers  shown  in  Fig  8(b).  Assuming  that 
the  dominant  TEM-mode  is  obliquely  incident,  the  electric  field  in  waveguide  A  can  be  written  as 

Eff  =  Aq  exp[-jkQ(z cos(^)  -I-  xsin(^))] ,  (11a) 


where  Aq 
incidence 
equations 


is  the  amplitude  of  the  incident  wave;  is  the  free  space  wave  number,  and  (j)  is  the  oblique 
angle.  The  magnetic  field  components  of  the  incident  TEM-mode  are  defined  using  Maxwell's 
as 


r  1 

jko\lQ  dz 


(lib) 


and 


1  ^Ey 

jkQp.Q  dx 


(11c) 


where  is  the  free  space  permeability.  Based  on  Eqs.  (11a)  and  (11b),  the  x-component  of  the  incident 
wave’s  magnetic  field  can  be  expressed  as 


=-Ao^^^^exp[-;A:o(zcos((|))-)-xsin(<t)))]  . 
M-O 


(12) 


Reflection  occurs  at  the  interface  between  waveguide  A  and  the  matching  network.  The  angle  of 
reflection  is  equal  to  the  angle  of  incidence.  The  reflected  wave’s  field  components  in  waveguide  A  are 

=  Rq  exp[ cos((|))  -  X  sin((l)))]  ( 1 3a) 

and 

=  /?o^^^^®^p[./^o(zcos((l))-xsin((|)))]  , 

M-O 


(13b) 
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where  7?q  is  the  amplitude  of  the  reflected  wave. 

Since  the  input  and  output  matching  networks  are  identical,  the  angle  at  which  the  plane  wave 
propagates  in  waveguide  C  (area  X  in  Fig.  8(b))  is  equal  to  the  angle  at  which  the  plane  wave  was 
incident  in  waveguide  A.  The  electric  and  magnetic  fields  of  the  dominant  mode  in  area  X  can  be  written 
as 

£(10)  _  exp[-;lo (zcos((t))  +  xsin((t)))]  (14a) 

and 

£r(i0)  _  ^^^exp[-7*o(^cos(<l))  +  ^sin((t)))]  ,  (14b) 

I^O 

where  is  the  amplitude  of  the  wave  in  waveguide  C  (area  X). 

In  area  IX  of  Fig.  8(b),  the  electromagnetic  fields  can  be  represented  as  a  sum  of  two  waves:  the  first 
wave  travels  in  the  positive  ^-direction,  with  an  amplitude  Tg;  and  the  second  wave  travels  in  the  negative 
z-direction,  with  an  amplitude  Rg.  The  fields  can  be  written  as 

E^y^=Tg  exp^-jf'A:^  (zcos((j))  +  xsin((|)))|  +  i?9exp[“jA:o(“ZCOs((|))  +  ;rsin((|)))]  (15a) 

and 

exp[“7A:o(zcos(<t))  +  xsin((|)))]+  /?9exp[--7/:o(”^^os((|))  +  :rsin((t)))]}  .  (15b) 

M'O 

Tg  and  Rg  are  to  be  determined  from  the  boundary  conditions  at  the  interfaces  between  areas  VIE  and  IX 
and  areas  IX  and  X.  The  oblique  angles  of  the  waves  propagating  in  area  IX  are  the  same  as  in  area  X 
because  both  areas  are  filled  with  the  same  dielectric  (air). 

The  boundary  conditions  require  that  the  tangential  components  of  electric  and  magnetic  fields  be 
continuous  at  all  discontinuity  interfaces.  The  interface  between  areas  IX  and  X  is  at 

9 

z  =  dg='^Li,  (16) 

1=1 

and  the  boundary  conditions  are 


(17a) 
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Substituting  electric  field  expressions  from  Eqs.  (14(a))  and  (15(a))  into  the  boundary  condition 
Eq.  (17(a))  followed  by  integration  over  the  interval  \—bQ/2<y<bQ/2\  leads  to 

Tio^o  exp[- jkQdg  cos((t))]  =  bg  \Tg  exp[-jkQdg  cos((l))]  +  Rg  exp[jkQdg  cos(<|))]}  .  (18a) 

Similarly,  substituting  magnetic  field  expressions  from  Eqs.  (14(b))  and  (15(b))  into  the  boundary 
condition  Eq.  (17(b))  followed  by  integration  over  the  interval  [-bg /2  <  y  <  ^>9/2]  leads  to 

-Tiq  ex]p[-jkQdg  cos(<|))]  =  -Tg  exp[-7^o^9  cos((|))]  +  Rg  exp[y^o^9  •  (18b) 

Dividing  Eq.  (18(a))  by  Eq.  (18(b)),  we  get 

^9  =  TgQg  exp[-j^o2^9  cos((|))]  ,  (19) 

where 


Using  Eqs.  (19)  and  (20),  the  electromagnetic  field  components  and  can  be  expressed  as 
=  Tg  exp[-jko (x sin((|))  +  dg  cos((t)))]  {exp[-;7:o (z-dg) cos((t))]  +  Qg  exp[#o (z  - ^^9 )cos((l))]} (21a) 
and 

=  79  ^^^^exp[-y7:o(xsin(<|))  +  ^9  cos((j)))] 

^0 

•{-exp[-7A:o(z-£?9)cos(<t))]  +  g9exp[y*o(^“^9)cos(<t))]}  •  (21b) 

In  area  VIII  of  Fig.  8(b)  the  field  components  can  be  written  as 

=  Tg  exp[-7A:o(^cos(<l))  +  xsin(<j)))]  +  Eg  exp[-jA:o(-^cos(<t))  +  xsin((j)))]  (22a) 

and 

^(g)  ^  cosW  r_j  sin(<t)))]  +Rg  exp  (“7^0  (“^  cos((|))  +  x  sin((t)))]}  .  (22b) 

Fo 
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The  first  term  corresponds  to  the  wave  traveling  in  the  positive  z-direction,  with  an  amplitude  Tg,  and  the 
second  term  corresponds  to  the  wave  traveling  in  the  negative  z-direction,  with  an  amplitude  /?g.  Tg  and 
/?g  are  to  be  determined  from  the  boundary  conditions  for  tangential  field  components  at  the  interface 
between  areas  Vni  and  IX: 

8 

z  =  J8^X^<-  (23) 

«=l 

The  boundary  conditions  are 

^<yA 

2  2 

2~^  2  (24a) 

^<y<-^ 

2  2 


and 

—  <y<—  .  (24b) 

2  2 

Substituting  electric  field  expressions  from  Eqs.  (21(a))  and  (22(a))  into  the  boundary  conditions  of 
Eq.  (24(a))  and  integrating  over  the  interval  [-bg  /2<y<bg  /2] ,  we  arrive  at 

Tg  expi-jkQdg  cos(<|)))»{exp(7A:oL9  cos((j)))  +  (29  exp(-y^0^9  cos((l)))}fe9 

=  [Tg  cxp(-jkod^  cos((|)))  +  /?g  exp(;^o^^8  cos(<t)))]fcg  .  (25a) 

Similarly,  substituting  magnetic  field  expressions  from  Eqs.  (21(b))  and  (22(b))  into  the  boundary 
condition  Eq.  (24(b))  and  integrating  over  the  interval  [-feg  /2  <  y  <  ftg  /2] ,  we  arrive  at 

T9  exp(-jkQdg  cos((|)))  •  [-exp(7^oL9  cos((|)))  -I-  Qg  expi-jkoLg  cos((t)))] 

=  — Tg  exp(-7^Qrfg  cos((|)))  + /?g  expljkgd^  cos(<t>))  .  (25b) 


f(9)_ 

Jt-y 


)  - 

0 


Dividing  Eq.  (25(a))  by  (25(b))  we  arrive  at 

^8  =  ^8*28  exp(-7A:o2Jg  cos((t)))  ,  (26a) 

where 


fg  + 


h 


<28  = 


f9-y 

bg 


1  +  09  ^^Pi-jkQ2Lg  cos((|))) 
-1  +  Qg  exp(- jkQ  2Lg  cos((t))) 


(26b) 


(26c) 
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Substitute  Eq.  (26(a))  into  Eqs.  (22(a))  and  (22(b)),  then 

exp[-jkoixsm(^)  +  d^  cos((|)))]{exp[-;A:o "  ^^8 ) cos((j))] 

+Qg  exp[jkQ(z-dg)cos(^)]} 

and 

^(8)  _  cos^ exp[-y*:o (x sin((|))  +  dg  cos((|)))]{- exp[-jko (z-dg) cos(<|))] 
M-O 

+08  ey^p[jko{z-dg)cosi(^)]}  . 

In  area  VH  in  Fig.  8(b),  the  field  components  can  be  written  as 

=  T-j  exp[-;A:o  sin((|))  +  V7Z)]  +  R-j  exp  [-7^0  (•^  sin((|))  -  vyz)] 


(27a) 


(27b) 


(28a) 


where 


hO)  ^  ZL exp[- jkQ  (x sin((|))  +  vyz)]  +  Rj  exp[-jkQ  (x sin((t))  -  V7z)]}  , 
^0 


V7 


^jej  -sin^((|))  . 


(28b) 

(28c) 


Continuing  this  process  of  writing  the  transverse  electromagnetic  field  components  in  each  consecutive 
area  and  applying  boundary  conditions  at  each  interface,  we  arrive  at  the  expression  for  the  dominant 
mode  reflection  coefficient  in  waveguide  A: 


where 


Rn  = 


h 


h 


Qi  = 


h 

\ 


Vm) 


Ui'+l  + 


7  ^ 


^•+1 


y^i+\ ) 


fi  = 


l  +  e,-  exp[-7fcov,-2L,  ] 

-l  +  Q;  exp[-7/!:oV/2L,  ]  ’ 


n  -sin^(<t))  ,  and 

/m  =-l>and 
=cos((l)). 


(29) 


(30) 


(31) 

(32) 

(33) 

(34) 

(35) 
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where  M  =  10.  To  determine  the  reflection  coefficient  for  the  matching  network  which  uses  four 
transformers,/!  is  needed.  It  is  computed  iteratively  as  follows.  Initially, /iq  =  -1  is  used  to  compute  Qg, 
and  then/9  is  computed;  fg  is  used  to  compute  Qg,  and  so  on  until/!  can  be  computed. 

A  similar  solution  was  derived  for  the  matching  network,  which  uses  three  transformers  shown  in 
Fig.  8(a).  It  is  expressed  by  the  same  set  of  equations  as  Eqs.  (29)-(35)  except  now  with  M=  8. 

Optimization  and  Numerical  Results 

Using  the  analytical  solutions  that  were  developed  in  the  previous  section,  optimization  algorithms 
based  on  the  minimax  approach  have  been  devised  that  minimize  the  maximum  value  of  the  reflection 
coefficient.  Computer  programs  (HOPT12,  HINC12,  HOPT22,  and  HINC22)  have  been  written  to 
implement  these  algorithms;  a  listing  of  these  programs  appears  in  the  Appendices.  The  programs 
optimize  the  length  of  each  transformer.  For  the  transformers  that  use  empty  parallel  plate  sections,  the 
height  is  also  optimized.  For  the  last  two  transformers  that  use  dielectric  materials,  the  programs  also 
optimize  their  dielectric  constants;  we  use  only  the  dielectric  constants  of  materials  that  can  readily  be 
obtained  from  manufacturers.  The  height  of  the  last  two  transformers  is  the  same  as  the  height  of  the 
guide  loaded  with  the  ferroelectric  material. 

The  programs  provide  the  optimum  values  of  the  matching  network  parameters  that  minimize  the 
magnitude  of  the  reflected  field  for  a  range  of  incidence  angles  over  the  frequency  band  of  interest  and  for 
the  required  variation  of  the  ferroelectric  dielectric  constant  to  obtain  360°  differential  phase  shift. 
HOPT12  provides  the  parameters  of  the  optimum  transformers  for  the  three-transformer  matching 
network  of  Fig.  8(a);  these  parameters  are  used  by  HINC12  to  compute  the  reflection  coefficient  at  each 
frequency,  incidence  angle,  and  dielectric  constant  of  interest.  HOPT22  and  HINC22  perform  similar 
computations  for  the  four-transformer  matching  network  of  Fig  8(b). 

The  computer  programs  HOPT12  and  HINC12  were  used  to  design  a  matching  network  with  three 
transformers  (Fig.  8(a)).  Three  scan  angle  ranges  were  considered:  0°  to  50°,  0°  to  55°,  and  0°  to  60°.  For 
each  scan  angle  range,  the  optimization  was  performed  over  two  different  frequency  bandwidths  centered 
at  10  GHz:  20%  (9  to  1 1  GHz)  and  30%  (8.5  to  1 1.5  GHz).  It  was  assumed  that  bQ  =  0.4  in.,  b2=b^  =  b^ 
=  0.1  in.,  and  E!  =£().  The  ferroelectric  dielectric  constant  was  assumed  to  be  100,  and  its  variation  was 
not  used  for  optimization.  The  minimax  optimization  minimizes  the  reflection  coefficient  over  the 
specified  frequency  range  and  scan  angles.  Table  2  presents  the  results  for  each  case.  For  each  case,  the 
reflection  coefficient  will  be  no  greater  than  the  value  provided  in  the  table  for  the  specified  scan  angle 
and  frequency  range.  If  the  optimization  was  performed  over  a  range  of  ferroelectric  dielectric  constants, 
then  the  reflection  coefficients  would  be  slightly  larger. 


Table  2  —  The  Maximum  Reflection  Coefficient 
for  a  Three-Transformer  Matching  Network 


Scan  Angle 

20%  Bandwidth 
(9  to  11  GHz) 

30%  Bandwidth 
(8.5  to  11.5  GHz) 

0°  to  50° 

0.240 

0.315 

0°  to  55° 

0.300 

0.396 

0°  to  60° 

0.400 

0.504 
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It  is  evident  from  Table  2  that  the  matching  network  that  uses  three  transformers  does  not  provide  a 
low  enough  reflection  coefficient  for  a  wide  range  of  incidence  angles.  However,  for  normal  incidence, 
the  reflection  coefficient  is  sufficiently  low.  The  optimization  for  normal  incidence  was  conducted 
assuming  that  the  dielectric  constant  of  the  ferroelectric  varies  from  120  to  80  (33%  tunability)  over 
a  frequency  range  of  8  to  12  GHz  (40%  bandwidth),  and  that  bQ  =  0.4  in.,  b2  =  b'^  =  =  0.1  in.,  and 

Ej  =  Eq.  Figure  9  shows  the  computed  results  for  normal  incidence.  The  corresponding  optimum 
parameters  are:  Lj  =  0.2956  in.,  L2  =  0.1860  in.,  1.3  =  0.0505  in.,  bi  =  0.2345  in.,  E2=  2.54,  and  E3=  35. 
These  results  were  compared  with  the  results  obtained  by  using  commercially  available  software 
“Touchstone”  for  a  rectangular  waveguide.  When  the  width  of  the  rectangular  waveguide  is  chosen  to  be 
greater  than  4^q,  the  results  obtained  by  Touchstone  are  identical  to  results  obtained  from  our  computer 
program  for  the  parallel  plate  waveguide  structure  under  consideration. 

1  - 
0.9 
0.8  - 

^  0.7  - 

.2 

£  0.6  - 

(L> 


Fig.  9  —  Reflection  coefficient  for  normal  incidence  with  three-section  matching  network 
for  various  dielectric  constants  (e^  =  80,  90,  100,  110,  and  120) 


The  matching  network  that  uses  four  transformers  (Fig.  8(b))  provides  sufficiently  low  reflection 
coefficient  even  for  the  case  of  oblique  incidence.  The  optimization  was  performed  using  the  computer 
programs  HOPT22  and  HINC22  for  two  operating  frequency  bands;  9.0  to  11.0  GHz  (20%  bandwidth) 
and  8.5  to  11.5  GHz  (30%  bandwidth).  Again,  three  scan  angle  ranges  were  considered:  0°  to  50°,  0°  to 
55°,  and  0°  to  60°,  and  it  was  assumed  that;  bQ  =  0.4  in.,  Z?3  =  I>4  =  =  0.1  in.,  and  Ej  =  E2  =  ^.  The 

unbiased  (maximum)  dielectric  constant  of  the  ferroelectric  was  assumed  to  be  100,  and  three  different 
tunabilities  were  used  in  the  optimization:  35%,  40%,  45%.  Knowing  that  ~  100,  for  each 

tunability,  E^  can  be  calculated  from  Eq.  (6).  E^  is  65,  60,  and  55  for  tunability  of  35%,  40%,  and 
45%,  respectively.  To  save  time,  optimization  was  performed  only  at  the  average  dielectric  constant  for 
each  tunability;  E^  is  82.5,  80,  and  77.5  for  tunability  of  35%,  40%,  and  45%,  respectively.  Variation 
of  the  dielectric  constant  was  not  used  in  the  optimization,  although  it  can  be  included.  The  results  are 
provided  in  Table  3. 
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Table  3  —  The  Maximum  Reflection  Coefficient  for  a  Four-Transformer  Matching  Network 


Scan  Angle 

20%  Bandwidth  (9  to  11  GHz) 

30%  Bandwidth  (8.5  to  1 1.5  GHz) 

35% 

Tunability 

40% 

Tunability 

45% 

Tunability 

35% 

Tunability 

40% 

Tunability 

45% 

Tunability 

0°  to  50° 

0.0700 

0.0558 

0.0538 

0.1190 

0.1035 

0.0960 

0°  to  55° 

0.0922 

0.0786 

0.0741 

0.1432 

0.1269 

0.1193 

0°  to  60° 

0.1366 

0.1144 

0.1227 

0.1918 

0.1769 

0.1808 

For  all  three  scan  angles,  the  case  with  35%  tunability  provides  the  largest  reflection  coefficient.  This 
is  because  the  variation  of  the  dielectric  constant  was  not  used  in  the  optimization.  Instead,  the  average 
was  used,  and  the  35%  tunability  case  has  the  largest  average  of  82.5.  As  it  was  explained  earlier,  the 
larger  the  e^,  the  more  difficult  the  matching  problem  becomes.  It  is  also  not  surprising  that  the  0°  to  60° 
scan  angle  case  provides  the  largest  reflection  coefficient  for  both  frequency  bandwidths.  If  the  variation 
of  would  have  been  used,  the  largest  reflection  coefficient  would  have  been  obtained  for  the  case  with 
30%  bandwidth,  0°  to  60°  scan  angle,  and  45%  tunability.  Next,  we  will  use  the  35%  tunability  and  0°  to 
60°  scan  angle  cases  to  show  that  variation  of  e^.  has  relatively  small  effect  on  the  reflection  coefficient. 

For  the  case  with  35%  tunability  i£r,avg  —  82.5),  0°  to  60°  scan  angle  and  20%  bandwidth,  the 
optimization  is  performed  for  a  dielectric  constant  of  82.5.  The  optimum  parameters  are:  Lj  =  0.3736  in., 
^2=  0.4186  in.,  L3  =  0.2116  in.,  L4  =  0.0539  in.,  =  0.2913  in.,  *2  =  0-1627  in.,  83  =  2.54,  and^=  30. 
Figure  10(a)  shows  that  the  reflection  coefficient  is  less  than  0.1366,  and  seven  curves  are  plotted 
corresponding  to  seven  incidence  angles.  With  all  other  parameters  kept  constant,  the  dielectric  constant 
of  the  ferroelectric  was  changed  to  65  and  100  to  examine  the  effect  of  applying  DC  bias  voltage  on  the 
reflection  coefficient.  Figures  10(b)  and  10(c)  show  the  results;  there  is  no  significant  change  in  the 
magnitude  of  the  reflection  coefficient  compared  to  that  of  Fig.  10(a). 

For  the  same  case,  when  the  optimization  is  performed  for  a  dielectric  constant  of  82.5  over  30% 
bandwidth,  the  optimum  parameters  are:  Lj  =  0.3685  in.,  L2  =  0.4048  in.,  L3  =  0.2074  in.;  L4  =  0.0531  in., 
Z>1  =  0.2991  in.,  ^2  =  0-1631  in.,  83  =  2.54,  and  8^  =  30.  Figure  11(a)  shows  that  the  reflection  coefficient  is 
less  than  0.1918  for  incidence  angles  up  to  60°.  As  before,  with  all  other  parameters  kept  constant,  the 
dielectric  constant  of  the  ferroelectric  was  changed  to  65  and  100  to  simulate  the  effect  of  applying  DC 
bias  voltage.  Figures  11(b)  and  11(c)  show  the  results.  Again,  there  is  no  significant  change  in  the 
magnitude  of  the  reflection  coefficient  compared  to  that  of  Fig  1 1(a).  These  results  indicated  that  the  lens 
is  well  matched  even  after  applying  DC  bias  voltage  to  scan  the  beam. 

Theoretical  results  presented  in  this  section  indicate  that  a  ferroelectric  lens  can  be  designed  that  is 
well  matched  over  large  scan  angles  and  broad  frequency  range.  The  reflection  coefficient  can  be  reduced 
further  (or  the  matching  can  be  improved)  by  decreasing  the  frequency  bandwidth  and/or  the  range  of 
incidence  angles  over  which  the  lens  operates.  If  necessary,  the  matching  can  also  be  improved  by 
increasing  the  number  of  transformers  in  the  matching  network. 
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Fig.  10(a)  —  Reflection  coefficient  for  several  incidence  angles  =  0,  10,  20,  30, 40,  50,  60°), 
35%  tunabUity  and  20%  bandwidth  with  four-section  matching  network 


Fig,  10(b)  —  Reflection  coefficient  for  several  incidence  angles  =  0,  10,  20,  30,  40,  50,  60°), 
35%  tunability  and  20%  bandwidth  with  four-section  matching  network 
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Fig.  10(c)  —  Reflection  coefficient  for  several  incidence  angles  ((j)  =  0, 10, 20,  30, 40,  50,  60°), 
35%  tunability  and  20%  bandwidth  with  four-section  matching  network 


EXPERIMENTAL  RESULTS 

Figure  12  shows  the  experimental  setup.  A  network  analyzer  is  connected  to  a  WR90  (X-band 
rectangular)  waveguide  through  a  coax-to-waveguide  transition.  The  width  of  the  WR90  is  flared  from 
0.9  in.  to  5  in.  to  create  an  fl^-plane  horn,  which  feeds  the  parallel  plate  column.  The  column  uses 
the  matching  network  with  three  transformers.  The  transformer  parameters  were  given  in  the  previous 
section,  and  Fig.  9  shows  the  reflection  coefficient.  The  network  analyzer  measures  the  .s-parameters 
of  the  column  as  a  function  of  frequency  and  DC  bias  voltage.  The  ferroelectrics  were  1  in.  long 
(in  the  direction  of  propagation),  0.05  in.  high,  and  5  in.  wide.  Experiments  were  performed  with  two 
different  ferroelectric  compositions,  Bao  ^QSrQ  4oTi03  with  60%  oxide  (abbreviated  60/60)  and 
Buq  55SrQ  45Ti03  with  60%  oxide  (abbreviated  55/60).  These  two  materials  offered  a  good  compromise 
among  8^,  tan  6,  and  tunability.  Table  1  shows  properties  of  these  materials;  =  100  and  tan  6  =  0.0079 
for  55/60;  and  =  118  and  tan  5  =  0.0129  for  60/60.  Figure  13  shows  the  measured  transmission  and 
reflection  coefficients  of  these  two  materials.  Since  the  input  and  output  reflection  coefficients  are  almost 
identical,  only  one  is  shown.  The  reflection  coefficient  is  sufficiently  small  over  a  wide  frequency  band  as 
the  theory  had  predicted  in  Fig.  9.  Because  of  the  higher  8^  and  the  higher  tan  6,  60/60  composition  has  a 
higher  loss. 

Figure  13  also  shows  that  the  loss  increases  with  frequency.  This  is  due  to  two  reasons;  first,  tan  8 
increases  with  frequency,  which  is  expected  for  ceramics;  second,  the  electrical  length  (in  terms  of 
wavelengths)  of  the  ferroelectric  in  the  direction  of  propagation  increases  with  frequency  since  the 
physical  length  is  kept  constant  (1  in.). 
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Fig.  11(a)  —  Reflection  coefficient  for  several  incidence  angles  (0  =  0, 10, 20, 30, 40, 50,  60°), 
35%  tunability  and  30%  bandwidth  with  four-section  matching  network 
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Fig.  1 1(b)  —  Reflection  coefficient  for  several  incidence  angles  ((()  =  0,  10,  20,  30, 40,  50,  60°), 
35%  tunability  and  30%  bandwidth  with  four-section  matching  network 
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Fig.  11(c)  — Reflection  coefficient  for  several  incidence  angles  =  0, 10, 20,  30, 40, 50,  60°), 
35%  tunability  and  30%  bandwidth  with  four-section  matching  network 


Fig.  12  —  Experimental  setup  for  one  column 
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As  discussed  earlier,  to  minimize  the  dielectric  loss  for  obtaining  360°  differential  phase  shift,  one 
would  like  to  obtain  as  high  a  tunability  as  possible.  Since  tunability  is  a  linear  function  of  the  DC  electric 
field  in  the  ferroelectric,  it  is  desirable  to  apply  as  high  a  DC  electric  field  as  possible.  If  the  dimension 
across  which  the  DC  bias  voltage  is  applied  is  kept  constant,  then  as  the  DC  bias  voltage  is  increased,  the 
DC  electric  field  increases.  Therefore,  during  experiments,  bias  voltage  was  increased  gradually  until 
dielectric  breakdown  occurred.  Figures  14(a)  and  14(b)  show  the  measured  phase  shift  of  the  55/60  and 
60/60  ferroelectrics  as  a  function  of  the  bias  voltage  for  various  frequencies.  Because  of  the  higher 
barium  content,  the  tunability  of  the  60/60  material  is  higher,  so  the  measured  phase  shift  is  higher.  As 
expected,  the  phase  shift  increases  linearly  with  frequency  because  the  electrical  length  of  the  material 
increases  with  frequency.  Since  ferroelectrics  are  good  insulators,  the  DC  current  requirements  are  very 
low.  For  example,  at  10  kV  bias  voltage,  the  DC  current  drawn  was  0.05  mA,  so  the  DC  power  dissipated 
is  only  0.5  W.  The  bias  voltage  can  be  reduced  by  further  bifurcating  the  ferroelectrics  (using  interdigital 
electrodes)  as  shown  in  Fig.  15.  Bifurcation  reduces  the  dimension  across  which  the  bias  voltage  is 
applied;  consequently,  the  required  bias  voltage  is  reduced  to  provide  the  same  DC  electric  field. 

As  DC  bias  voltage  is  applied  to  the  ferroelectrics,  the  dielectric  constant  decreases,  which  is  the 
basis  for  obtaining  the  differential  phase  shift,  but  it  also  complicates  the  impedance  matching  problem. 
However,  our  design  optimization  is  able  to  cope  with  this  difficult  problem.  Figures  16(a)  and  16(b) 
show  the  reflection  coefficients  for  the  55/60  and  60/60  materials  as  a  function  of  frequency  for  various 
bias  voltages.  The  same  parallel  plate  waveguide,  with  the  same  dielectric  transformers,  is  used  for  both 
materials.  The  standing  wave  ratio  (SWR)  is  less  than  2  for  frequency  range  of  8  to  12  GHz,  as 
theoretically  predicted  earlier. 
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Pig^  15  —  Scheme  to  reduce  the  bias  voltage 


Considerable  experimentation  was  performed  to  gain  experience  with  high-dielectric-constant  loaded 
waveguide.  Air  gaps  between  the  ferroelectric  and  the  parallel  conducting  plates  cause  considerable 
problems  with  resonances.  Poor  electrodes  will  cause  the  electric  field  to  concentrate  in  the  air  gaps 
causing  breakdown;  also,  the  field  in  the  ferroelectric  is  reduced,  which  decreases  the  ferroelecmc 
tunability.  The  ferroelectrics  are  supplied  by  ARL  with  the  top  and  bottom  surfaces  metellized  with  a 
baked-on  silver  ink.  When  placed  in  the  parallel  plate  waveguide,  the  top  of  the  ferroelectnc  is  m  contact 
with  the  top  parallel  plate  at  ground  potential,  and  the  bottom  of  the  ferroelectnc  is  in  contact  wU  m 
electrode  that  is  connected  to  the  high-voltage  power  supply.  On  the  other  side  of  this  electrode  is  the 
other  piece  of  the  ferroelectric  whose  top  is  in  contact  with  the  high-voltage  electrode  and  whose  bottom 
is  in  contact  with  the  bottom  parallel  plate  at  ground  potential.  The  small  air  gaps  between  the 
ferroelectric  and  the  parallel  plates  are  filled  using  thin  indium  sheets  and  various  types  of  conductive 
greases  and  caulks.  Indium  is  used  because  it  has  some  ductility.  It  was  discovered  that  air  gaps  as  small 
as  a  fraction  of  a  mil  can  cause  problems. 


The  experimental  results  presented  here  agree  well  with  the  theoretical  results  of  the  lens  aperture 
matching.  Control  of  the  column  phase  shift  with  the  applied  DC  voltage  was  also  demonstrated 
experimentally.  In  summary,  the  experimental  results  on  one  column  of  the  lens  proved  that  t  e 
ferroelectric  lens  concept  is  viable  and  sound.  Work  is  continuing  in  building  and  demonstrating  a  sma 

lens. 


6.  CONCLUSIONS 

A  new  concept  for  low-cost  phased  arrays  using  a  voltage  controlled  ferroelectric  lens  was  proposed. 
The  ferroelectric  lens  phased  array  uniquely  incorporates  bulk  phase  shifting  (the  array  does  not  contain 
discrete  phase  shifters)  using  ferroelectric  material,  which  reduces  significantly  the  number  of  phase 
shifters  and  drivers.  Furthermore,  row-column  beam  steering  is  used.  These  features,  m  turn,  reduce  the 
cost  of  the  phased  array.  This  report  presented  the  ferroelectric  lens  concept,  its  theoretical  analysis, 
design  and  experimental  results.  The  theoretical  results  agree  well  with  experimental  results,  proving  that 
the  ferroelectric  lens  concept  is  very  practical.  It  was  shown  that  the  lens  surface  can  be  matched  very 
well  to  free  space  over  a  broad  range  of  frequencies  and  large  scan  angles.  A  number  of  low-cost  phased 
array  configurations  using  the  ferroelectric  lens  were  also  described. 
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Appendix  A 

COMPUTER  PROGRAM  “HOPT12” 


PROGRAM  DESCRIPTION 

“HOPT12”  is  an  optimization  software  for  the  impedance  matching  network  of  Fig.  8(a),  which  uses 
three  transformers.  It  optimizes  the  length  of  all  three  transformers.  The  height  of  the  first  transformer  is 
also  optimized.  For  the  last  two  transformers,  the  dielectric  constants  are  also  optimized.  For  the  second 
transformer,  the  dielectric  constants  that  are  used  for  optimization  are  1,  2.54,  and  4.3  since  these  are 
dielectrics  that  we  can  readily  obtain.  Polystyrene  has  a  dielectric  constant  of  2.54,  and  it  has  always 
turned  out  to  be  the  optimum  material  for  the  second  transformer.  For  the  third  transformer,  we  use  the 
following  dielectric  constants:  11.6,  13.0,  14.0,  15.0,  16.0,  18.0,  20.0,  30.0,  35.0,  38.0,  40.0,  50.0,  55.0, 
and  70.0.  Again,  these  values  are  chosen  because  we  can  readily  purchase  ceramics  with  these  dielectric 
constants. 

The  optimization  is  performed  for  all  combinations  of  the  dielectric  constants  in  the  second  and  third 
transformers.  Input  data,  which  will  be  requested  by  the  program,  are 

a.  heights  of  the  waveguides  A,  B,  and  C; 

b.  length  of  the  waveguide  B  (length  of  the  ferroelectric);  waveguides  A  and  C  are  assumed  to  be 
semi-infinite; 

c.  maximum  (at  zero  bias)  and  minimum  (at  maximum  bias)  dielectric  constant  of  the  ferroelectric 
and  the  step  increment  to  be  used  for  optimization  as  the  dielectric  constant  is  changed  from  its 
maximum  to  minimum  value; 

d.  maximum  scan  angle  and  a  step  increment; 

e.  upper  and  lower  frequency  and  a  step  increment;  and 

f.  value  of  the  dielectric  constant  of  the  material  in  the  second  transformer. 

All  lengths  are  to  be  entered  in  inches,  frequencies  in  GHz,  and  angles  in  degrees. 

The  output  of  the  program  is  written  into  the  file  "bbbbb.dat".  The  output  consists  of  data  groups  for 
each  combination  of  the  second  and  third  transformer  dielectric  constants.  The  output  parameters  are 
denoted  as:  bl,  the  height  of  the  first  transformer;  and  LI,  L2  and  L3,  length  of  each  transformer.  For 
each  combination,  the  program  also  provides  the  upper  limit  for  the  reflection  coefficient,  which  is 
denoted  in  the  program's  output  as  "Fvalue"  and  listed  in  Table  2  for  the  cases  that  we  had  considered. 
The  matching  network  with  the  smallest  "Fvalue"  is  the  optimum  one. 
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LISTING  OF  PROGRAM 

c  3  transformers  optimization 


common  /blkl/  iangle,anglestp,frl,fru,frs,eps2,eps3,allls 

common  /blk2/  bhl,bhs 

common  /blk3/  epssin,epssend,epssstep 

integer  n 

parameter  (n=4) 

integer  k,maxfcn,nout,i502,i503,i504,i505,i506 
realaaa,hinciden,ftol,fvalue,s,x(n),xguess(n),ai507 

real  ep3(14),zero 

external  hincide3,umach,umpol 

open  (unit=l  ,file='ccccc.dat',status='new') 

print  *,  "Results  in  file:  ccccc.dat" 

zero=0.0 

ep3(l)=11.6 

ep3(2)=13.0 

ep3(3)=14.0 

ep3(4)=15.0 

ep3(5)=16.0 

ep3(6)=18.0 

ep3(7)=20.0 

ep3(8)=30.0 

ep3(9)=35.0 

ep3(10)=38.0 

ep3(ll)=40.0 

ep3(12)=50.0 

ep3(13)=55.0 

ep3(14)=70.0 

print  *,  "Enter  heights  of  large  and  small  waveguides  (inches)" 
read  *,  bhl,bhs 

print  *,  "Enter  length  of  smaller  waveguide  (inches)" 
read  *,  allls 

print  *,  "Enter  init.,  final,  and  step  of  dielectric  constant" 

print  *,  "of  smaller  waveguide" 

read  *,  epssin,epssend,epssstep 

print  *,  "Enter  maximium  oblique  scan  angle  and  step" 

read  *,  angle,anglestp 

iangle=angle/anglestp 

print  *,  "Enter  ffeq.  low  boundary,  upper  boundary,  freq.  step" 
read  *,  frl,frii,ffs 

print  *,  "Enter  dielectric  constant  of  second  transformer" 

read  *,  eps2 

do502i502=l,14 

fv=1000.0 

eps3=ep3(i502) 

do  503  i503=l,3 

xguess(2)=0.1  +0. 1  *(i503- 1 ) 

do  504  i504=l,3 


Voltage-Controlled  Ferroelectric  Lens  Phased  Array 


33 


xguess(3)=0.1+0.1*(i504-l) 
do  505  i505=l,2 

xguess(4)=0.045+0.025*(i505-l) 

ai506=0.05 

ii5 06=INT((0.4-0. 1  )/ai506)-4 

do  506  i506=l,ii506 

xguess(l)=0.15+ai506*(i506-l) 

ftol=1.0e-4 

maxfcn=3000 

s=1.0 

call  umpol  (hmcide3,n,xguess,s,ftol,maxfcn,x,fvalue) 

if  (fvalue.LT.fv.AND.x(2).GT.zero.AND.x(3).GT.zero.AND. 

$x(4).GT.zero.AND.x(4).LE.0.08 

$.AND.x(l).GE.0.1.AND.x(l).LE.0.4)then 

fv=fvalue 

bl=x(l) 

all=x(2) 

al2=x(3) 

al3=x(4) 

end  if 

506  continue 

505  continue 

504  continue 

503  continue 
write(l,778) 

778  formate  lx,' - ') 

write(l,777)  eps2,eps3,bl,all,al2,al3,fv 

777  fonnat(lx,'  eps2=',fl0.4/lx,'  eps3=',fl0.4/lx, 

$'  bl=',fl0.4/lx,'  Ll=',fl  0.4/lx,'  L2=',fl0.4/ 

$lx,'  L3=',fl0.4/lx,'  —  Fvalue — =',fl0.4) 

502  continue 
continue 
stop 
end 

subroutine  hincide3  (n,x,f) 
implicit  real*4  (a-h,o-z) 

common  ^lkl/  iangle,anglestp,frl,fru,frs,eps2,eps3,alllsl 

common  /blk2/  bhl  1  ,bhs  1 

common  /blk3/  epssinl,epssendl,epssstpl 

integer  n 

complex  qq,ff,ffteml,fftem,qqtem,ref 
dimension  x(n) 

dimension  qq(  1 00),ff(  1 00),bh(  1 00),al(  1 00) 

dimension  eps(100),anu(100),gam(l  00) 

one=1.0 

two=2.0 

ten=10.0 

api=acos(-one) 

api2=api/two 

apit=api 
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api2t=api2 
constl=0. 16933333 
valuel=0.0 
i234=iangle 

aepss=(epssend  1  -epssin  1  )/epssstp  1 

iepss=NINT  (aepss) 

do  30 130=1, iepss+1 

epss=epssin  1+epssstp  1  *(130- 1 ) 

do  20 120=1,1234+1 

phldeg=anglestp*(120- 1 ) 

phl=phldeg*  api/1 80.0 

bhl=bhll 

bhs=bhsl 

allls=alllsl 

nst=3 

nstl=2*nst+l 

nst2=2*nst+2 

bh(nst+l)=bhs 

eps(nst+l)=epss 

al(nst+l)=allls 

bh(l)=x(l) 

bh(2)=bhs 

bh(3)=bhs 

bh(5)=bh(3) 

bh(6)=bh(2) 

bh(7)=bh(l) 

al(l)=x(2) 

al(2)=x(3) 

al(3)=x(4) 

al(5)=al(3) 

al(6)=al(2) 

al(7)=al(l) 

eps(l)=1.0 

eps(2)=eps2 

eps(3)=eps3 

eps(5)=eps(3) 

eps(6)=eps(2) 

eps(7)=eps(l) 

frequ=fru 

freql=frl 

frstep=frs 

numfr=(frequ-freql)/frstep 

do  10  lfT=0,nuinfr 

fr=freql+frstep*lfr 

wnum=apl*constl  *fr 

cccc  1  =(bhybh(7)-one)/(bhl/bh(7)+one) 

qq(7)=cmplx(cccc  1,0.0) 

ffteml=cmplx(0.0,-wnum*two*al(7)*cos(phl)) 

fftem=exp(fftenil) 

ff(7)=(one+qq(7)*fftem)/(-one+qq(7)*fftem) 
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anu(6)=sqrt(eps(6)-(sin(phi))**2) 

qqtem=(bh(7)/bh(6))*(anu(6)/cos(phi))*ff(7) 

qq(6)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0,-wnuni*two*al(6)*anu(6)) 

fftem=exp(ffteml) 

ff(6)=(one+qq(6)*fftem)/(-one+qq(6)*fftem) 

anu(5)=sqit(eps(5)-(sin(phi))**2) 

qqtem=(anu(5)/anu(6))*ff(6) 

qq(5)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0,-wnuni*two*al(5)*anu(5)) 

fftem=exp(ffteml ) 

ff(5)=(one+qq(5)*fftem)/(-one+qq(5)*fftem) 

anu(4)=sqrt(eps(4)-(sin(phi))**2) 

qqtem=(anu(4)/anu(5))*ff(5) 

qq(4)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0,-wnum*two*al(4)*anu(4)) 

fftem=exp(ffteml ) 

ff(4)=(one+qq(4)*fftem)/(-one+qq(4)*fftem) 

anu(3)=sqrt(eps(3)-(siii(phi))**2) 

qqtem=(anu(3)/anu(4))*ff(4) 

qq(3)=(qqtem+one)/(qqtem-one) 

ffteml =cmplx(0.0,- wnum*two*al(3)  *anu(3)) 

fftem=exp(ffteml ) 

ff(3)=(one+qq(3)*fftem)/(-one+qq(3)*fftem) 

anu(2)=sqrt(eps(2)-(sm(phi))  *  *2) 

qqtem=(anu(2)/anu(3))*ff(3) 

qq(2)=(qqtem+one)/(qqtem-one) 

fifteml=cmplx(0.0,-wnum*two*al(2)*anu(2)) 

fftem=exp(ffteml) 

ff(2)=(one+qq(2)*ffl;eni)/(-one+qq(2)*fftem) 
qqtem=(cos(phi)/anu(2))  *ff(2) 
qq(l)=(qqtem+bh(l)/bh(2))/(qqtem-bh(l)/bh(2)) 
fftem  1 =cmplx(0.0,-wnum*two*al(  1 )  *cos(phi)) 
fftem=exp(fft:eml ) 

ff(  1  )=(one+qq(  1 )  *fftem)/(-one+qq(  1 )  *fftem) 
ref=(ff(l)+bhm(l))/(ff(l)-bhM)h(l)) 
refr=real(ref) 
refi=aimag(ref) 
refabs=sqrt(refr*  *2+refi*  *2) 
if  (refabs.GT.valuel)  then 
value  l=refabs 
end  if 
10  continue 
20  continue 
30  continu 
f=valuel 
return 
end 


Appendix  B 

COMPUTER  PROGRAM  “HINC12” 


PROGRAM  DESCRIPTION 

“HINC12”  uses  the  optimum  impedance  matching  network  parameters  generated  by  “HOPT12”  for 
the  matching  network  of  Fig.  8(a),  which  uses  three  transformers.  “fflNC12”  uses  these  parameters  to 
compute  the  reflection  coefficient  at  each  frequency,  incidence  angle,  and  dielectric  constant  of  interest. 
The  program  writes  outputs  of  its  calculation  into  the  data  file  "aaaaa.dat".  The  output  data  file  consists  of 
four  columns:  the  first  column  lists  the  frequency  in  GHz,  the  second  column  gives  the  absolute  value  of 
the  reflection  coefficient,  the  third  column  is  the  absolute  value  of  the  reflection  coefficient  in  dB,  and  the 
fourth  column  is  the  phase  of  the  reflection  coefficient  in  radians.  Input  data  to  the  program  are 

a.  heights  of  the  waveguides  A,  B,  and  C; 

b.  length  of  the  waveguide  B  (length  of  the  ferroelectric);  waveguides  A  and  C  are  assumed  to  be 
semi-infinite; 

c.  dielectric  constant  of  the  ferroelectric; 

d.  incidence  (scan)  angle; 

e.  upper  and  lower  frequency  and  a  step  increment; 

f.  dielectric  constant  of  the  material  in  the  second  and  third  transformer; 

g.  lengths  of  all  three  transformers;  and 

h.  height  of  the  first  transformer. 

All  lengths  are  to  be  entered  in  inches,  frequencies  in  GHz,  and  angles  in  degrees. 

LISTING  OF  PROGRAM 
c  symmetrical  transition:  3  transformers 


implicit  double  precision  (a-h,o-z) 

complex  qq,ff,ffteml,fftem,qqtem,ref 

dimension  qq(l  00),ff(100),bh(  1 00),al(  1 00) 

dimension  eps(  1 00),anu(  1 00),gam(  1 00) 

open  (unit=l,file='aaaaa.dat',status-new') 

print  *,  "Results  are  in  file:  aaaaa.dat" 

one=1.0d0 

two=2.0d0 

ten=10.0d0 

api=acos(-one) 

api2=api/two 

apit=api 

api2t=api2 

constl=0.16933333d0 
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print  *,  "Enter  incidence  angle  in  degrees" 
read  (5,*)  phideg 
phi=phideg*api/l  SO.OdO 

print  *,  "Enter  low  frequency,  upper  frequency,  freq.step" 
read  (5,*)  freql,frequ,ffstep 

print  *,  "Enter  larger  waveguide  height  (in  inches)" 
read  (5,*)  bhl 

print  *,  "Enter  smaller  waveguide  height  (in  inches)" 
read  (5,*)  bhs 

print  *,  "Enter  smaller  waveguide  length  (in  inches)" 
read  (5,*)  allls 

print  *,"Enter  smaller  waveguide  dielectric  constant  (relative)" 

read  (5,*)  epss 

nst=3 

print  *,  "Enter  height  of  the  first  transformer  (in  inches)" 
read  (5,*)  bh(l) 

print  *,  "Enter  length  of  the  first  transformer  (in  inches)" 
read  (5,*)  al(l) 

print  *,  "Enter  length  of  the  second  transformer  (in  inches)" 
read  (5,*)  al(2) 

print  *,  "Enter  length  of  the  third  transformer  (in  inches)" 

read  (5,*)  al(3) 

nstl=2*nst+l 

nst2=2*nst+2 

bh(nst+l)=bhs 

eps(nst+l)=epss 

al(nst+l)=allls 

bh(2)=bhs 

bh(3)=bhs 

bh(5)=bh(3) 

bh(6)=bh(2) 

bh(7)=bh(l) 

al(5)=al(3) 

al(6)=al(2) 

al(7)=al(l) 

print  *,  "Enter  dielectric  constant  of  the  second  transformer" 
read  (5,*)  eps(2) 

print  *,  "Enter  dielectric  constant  of  the  third  transformer" 

read  (5,*)  eps(3) 

eps(l)=1.0d0 

eps(5)=eps(3) 

eps(6)=eps(2) 

eps(7)=eps(l) 

numfr=(frequ-freql)/frstep 

do  10  ifr=0,numfr 

fr=ffeql+ffstep*ifr 

wnum=api*constl  *6" 

ccccl=(bhl/bh(7)-one)/(bhlA»h(7)+one) 

qq(7)=cmplx(cccc  1 ,0.0d0) 

ffteml=cmplx(0.0d0,-wnum*two*al(7)*dcos(phi)) 
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fftem=exp(ffteml) 

ff(7)=(one+qq(7)*fftem)/(-one+qq(7)*fftem) 
anu(6)=sqrt(eps(6)-(sin(phi))  *  *2) 
qqtem=(bh(7)/bh(6))*(anu(6)/cos(phi))*ff(7) 
qq(6)=(qqtem+one)/(qqtem-one) 
ffteml=cmplx(0.0d0,-wnum*two*al(6)*anu(6)) 
fftem=exp(fftem  1 ) 

ff(6)=(one+qq(6)*fftem)/(-one+qq(6)*fftem) 

anu(5)=sqrt(eps(5)-(sin(phi))**2) 

qqtem=(anu(5)/anu(6))*ff(6) 

qq(5)=(qqtem+one)/(qqtem-one) 

ffieml=cmplx(O.OdO,-wnum*two*al(5)*anu(5)) 

fftein=exp(ffteml) 

ff(5)=(one+qq(5)*fftem)/(-one+qq(5)*ffiem) 

anu(4)=sqrt(eps(4)-(sin(phi))**2) 

qqtem=(anu(4)/anu(5))*ff(5) 

<l‘l(4)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0d0,-wnum*two*al(4)*anu(4)) 

fftein=exp(ffteml ) 

ff(4)=(one+qq(4)*fftem)/(-one+qq(4)*fftem) 

anu(3)=sqrt(eps(3)-(sin(phi))**2) 

qqtem=(anu(3)/anu(4))*ff(4) 

qq(3)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(O.OdO,-wnum*two*al(3)*anu(3)) 

fftem=exp(ffteml ) 

ff(3)=(one+qq(3)*fftem)/(-one+qq(3)*fftem) 

anu(2)=sqrt(eps(2)-(sin(phi))  *  *2) 

qqtem=(anu(2)/anu(3))*ff(3) 

qq(2)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0d0,-wnum*two*al(2)*anu(2)) 

fftem=exp(ffteml) 

ff(2)=(one+qq(2)*fftem)/(-one+qq(2)*fftem) 

qqtem=(cos(phi)/anu(2))*ff(2) 

qq(  1  )=(qqtem+bh(  1  )/bh(2))/(qqtem-bh(  1  )/bh(2)) 

ffteml=cmplx(O.OdO,-wnum*two*al(l)*cos(phi)) 

fftem=exp(ffteml ) 

ff(l)=(one+qq(l)*fftem)/(-one+qq(l)*fftem) 
ref=(ff(l)+bhl^h(l))/(ff(l)-bhl^h(l)) 
refr=real(ref) 
refi=aimag(ref) 
refabs=sqrt(refr*  *2+refi*  *2) 
refdb=20.0d0*log  1 0(refabs) 
refarg=atan(refi/refr) 
write  (1,20)  fT,refabs,refdb,refarg 
20  fonnat(lx,f5.2,lx,f8.5,lx,f6.2, 

$lx,f8.3) 

10  continue 
stop 
end 


Appendix  C 

COMPUTER  PROGRAM  “HOPT22” 

PROGRAM  DESCRIPTION 

“HOPT22”  is  an  optimization  software  for  the  impedance  matching  network  of  Fig.  8(b),  which  uses 
four  transformers.  It  optimizes  the  length  of  all  four  transformers.  The  height  of  the  first  two  transformers 
is  also  optimized.  For  the  last  two  transformers,  the  dielectric  constants  are  also  optimized.  For  the  third 
transformer,  the  dielectric  constants  that  are  used  for  optimization  are  1,  2.54,  and  4.3  since  these  are 
dielectrics  that  we  can  readily  obtain.  Polystyrene  has  a  dielectric  constant  of  2.54,  and  it  has  always 
turned  out  to  be  the  optimum  material  for  the  third  transformer.  For  the  fourth  transformer,  we  use  the 
following  dielectric  constants:  11.6,  13.0,  14.0,  15.0,  16.0,  18.0,  20.0,  30.0,  35.0,  38.0,  40.0,  50.0,  55.0, 

and  70.0.  Again,  these  values  are  chosen  because  we  can  readily  purchase  ceramics  with  these  dielectric 
constants. 

The  optimization  is  performed  for  all  combinations  of  the  dielectric  constants  in  the  third  and  fourth 
transformers.  Input  data,  which  will  be  requested  by  the  program,  are 

a.  heights  of  the  waveguides  A,  B,  and  C; 

b.  length  of  the  waveguide  B  (length  of  the  ferroelectric);  waveguides  A  and  C  are  assumed  to  be 
semi-infinite; 

c.  maximum  (at  zero  bias)  and  minimum  (at  maximum  bias)  dielectric  constant  of  the  ferroelectric, 
and  the  step  increment  to  be  used  for  optimization  as  the  dielectric  constant  is  changed  from  its 
maximum  to  minimum  value; 

d.  maximum  scan  angle  and  a  step  increment; 

e.  upper  and  lower  frequency  and  a  step  increment;  and 

f.  value  of  the  dielectric  constant  of  the  material  in  the  third  transformer. 

All  lengths  are  to  be  entered  in  inches,  frequencies  in  GHz,  and  angles  in  degrees. 

The  output  of  the  program  is  written  into  the  file  "bbbbb.dat".  The  output  consists  of  data  groups  for 
each  combination  of  the  third  and  fourth  transformer  dielectric  constants.  The  output  parameters  are 
denoted  as:  bl  and  b2,  the  height  of  the  first  two  transformers,  and  LI,  L2,  L3  and  L4,  length  of  each 
transformer.  For  each  combination,  the  program  also  provides  the  upper  limit  for  the  reflection 
coefficient,  which  is  denoted  in  the  program's  output  as  "Fvalue",  and  listed  in  Table  3  for  the  cases  that 
we  had  considered.  The  matching  network  with  the  smallest  "Fvalue"  is  the  optimum  one. 
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LISTING  OF  PROGRAM 

c  4  transformers  optimization 


common  /blkl/  iangle,anglestp,frl,fru,frs,eps3,eps4,ams 

common  /blk2/  bhl,bhs 

common  /blk3/  epssin,epssend,epssstep 

integer  n 

parameter  (n=6) 

integer  k,maxfcn,nout,ii507,i502,i503,i504,i505,i506, 
$i507,i508 

realaaa,hinciden,ftol,fvalue,s,x(n),xguess(n),ai507 

real  ep4(14),zero 

external  hincide4,umach,umpol 

open  (unit=l  ,file='bbbbb.dat',status='new') 

print  *,  "Results  in  file;  bbbbb.dat" 

zero=0.0 

ep4(l)=11.6 

ep4(2)=13.0 

ep4(3)=14.0 

ep4(4)=15.0 

ep4(5)=16.0 

ep4(6)=18.0 

ep4(7)=20.0 

ep4(8)=30.0 

ep4(9)=35.0 

ep4(10)=38.0 

ep4(ll)=40.0 

ep4(12)=50.0 

ep4(13)=55.0 

ep4(14)=70.0 

print  *,  "Enter  heights  of  large  and  small  waveguides(inches)" 
read  *,  bhl,bhs 

print  *,  "Enter  length  of  smaller  waveguide  (inches)" 
read  *,  allls 

print  *,  "Enter  init,  final,  step  of  dielectric  constant" 
print  *,  "of  smaller  waveguide" 
read  *,  epssin,epssend,epssstep 

print  *,  "Enter  maximum  oblique  scan  angle  and  step" 

read  *,  angle, anglestp 

iangle=angle/anglestp 

print  *,  "Enter  freq.  low  boundary,  upper  boundary,  fi-eq.  step" 
read  *,  frl,fiTi,frs 

print  *,  "Enter  dielectric  constant  of  third  transformer" 

read  *,  eps3 

do  502  i502=l,14 

fv=l 000.0 

eps4=ep4(i502) 

do  503i503=l,3 

xguess(l  )=0. 1+0. 1  *(i503- 1 ) 
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do  504  i504=l,2 

xguess(2)=0.045+0.025 *  (i504- 1 ) 
do  505  i505=l,3 
xguess(3)=0.25+0.05*(i505-l) 
do  506  i506=l,2 
xguess(4)=0.3+0. 1  *(i506- 1) 
ai507=0.05 

ii507=INT((xguess(3)-0.1)/ai507)-2 

do  507  i507=l,ii507 

xguess(5)=0.  l+ai507*(i507- 1) 

do  508  i508=l,2 

xguess(6)=0.3+0.1*(i508-l) 

ftol=1.0e-4 

niaxfcn=3000 

s=1.0 

call  umpol  (hincide4,n,xguess,s,ftol,maxfcn,x,fvalue) 

if(fvalue.LT.fv.AND.x(3).GT.zero.AND.x(3).LT.bhl.AND. 

$x(5).GT.zero.AND. 

$x(4).GT.zero.AND.x(6).GT.zero.AND.x(l).GT.zero.AND. 

$x(2).GT.zero)  then 

fv=fvalue 

bl=x(3) 

b2=x(5) 

all=x(4) 

al2=x(6) 

al3=x(l) 

al4=x(2) 

end  if 

508  continue 
507  continue 
506  continue 
505  continue 
504  continue 
503  continue 
write(l,510) 

510fonnat(lx,' - ') 

write(l ,601)  eps3,eps4,bl ,b2,all ,al2,al3,al4,fv 
601  fomiat(lx,'  eps3=',fl0.4/lx,'  eps4=',f  10.4/lx, 

&'  bl=',fl0.4/lx,'  b2=',fl0.4/lx, 

&’  Ll=',fl0.4/lx,’  L2=’,fl0.4/ 

&lx,'  L3=',fl0.4/lx,'  L4=',fl0.4/lx, 

&' —  Fvalue  —  =',fl0.4) 

502  continue 
write(l,510) 
continue 
stop 
end 

subroutine  hincide4  (n,x,f) 
implicit  real*4  (a-h,o-z) 

common  /blkl/  iangle,anglestp,frl,fru,frs,eps3,eps4,alllsl 
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common  /blk2/  bhll,bhsl 

common  /blk3/  epssinl,epssendl,epssstpl 

integer  n 

complex  qq,ff,ffteml  ,fftem,qqtem,ref 
dimension  x(n) 

dimension  qq(  1 00),ff(  1 00)  ,bh(  1 00)  ,al(  1 00) 

dimension  eps(100),anu(100),gam(100) 

one=1.0 

two=2.0 

ten=10.0 

api=acos(-one) 

api2=api/two 

apit=api 

api2t=api2 

const 1  =0. 1 6933333 

value  1=0.0 

aepss=(epssendl  -epssinl  )/epssstp  1 

iepss=NINT(aepss) 

do  30  i30=l,iepss+l 

epss=epssinl+epssstpl*(i30-l) 

i234=iangle 

do  20 120=1,1234+1 

phideg=anglestp*(i20- 1) 

phi=phideg*api/l  80.0 

bhl=bhll 

bhs=bhsl 

allls=alllsl 

nst=4 

nstl=2*nst+l 

nst2=2*nst+2 

bh(nst+l)=bhs 

eps(nst+l  )=epss 

al(nst+l)=allls 

bh(l)=x(3) 

bh(2)=x(5) 

bh(3)=0.1 

bh(4)=0.1 

bh(6)=bh(4) 

bh(7)=bh(3) 

bh(8)=bh(2) 

bh(9)=bh(l) 

al(l)=x(4) 

al(2)=x(6) 

al(3)=x(l) 

al(4)=x(2) 

al(6)=al(4) 

al(7)=al(3) 

al(8)=al(2) 

al(9)=al(l) 

eps(l)=1.0 
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eps(2)=1.0 

eps(3)=eps3 

eps(4)=eps4 

eps(6)=eps(4) 

eps(7)=eps(3) 

eps(8)=eps(2) 

eps(9)=eps(l) 

frequ=fru 

fTeql=M 

frstep=frs 

numfr=(frequ-freql)/firstep 

do  10  ifT=0,numfr 

fi^freql+frstep*ifr 

wnuin=api*constl  *fr 

cccc  1 =(bhl/bh(9)-one)/(bhl/bh(9)+one) 

qq(9)=cmplx(cccc  1,0.0) 

ffteml=cmplx(0.0,-wnuni*two*al(9)*cos(phi)) 

fftem=exp(ffteml) 

ff(9)=(one+qq(9)*f¥tem)/(-one+qq(9)*fftem) 

anu(8)=sqrt(eps(8)-(sin(phi))**2) 

qqtein=(bh(9)/bh(8))*(anu(8)/cos(phi))*ff(9) 

qq(8)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0,-wnum*two*al(8)*anu(8)) 

fftem=exp(ffteinl) 

ff(8)=(one+qq(8)*fftem)/(-one+qq(8)*fftem) 
anu(7)=sqrt(eps(7)-(sin(phi))  *  *2) 
qqtem=(bh(8)/bh(7))*(anu(7)/anu(8))*ff(8) 
‘l<l(7)=(qqtem+one)/(qqtem-one) 
ffteml=cmplx(0.0,-wnum*two*al(7)*anu(7)) 
fftem=exp(fftem  1 ) 

ff(7)=(one+qq(7)*fftem)/(-one+qq(7)*fftem) 

anu(6)=sqrt(eps(6)-(sin(phi))**2) 

qqtem=(bh(7)/bh(6))*(anu(6)/anu(7))*ff(7) 

q9(6)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0,-wnum*two*al(6)*anu(6)) 

fftem=exp(ffteml) 

ff(6)=(one+qq(6)*fftein)/(-one+qq(6)*ffieni) 

anu(5)=sqrt(eps(5)-(sin(phi))**2) 

qqtem=(bh(6)/bh(5))*(anu(5)/anu(6))*ff(6) 

qq(5)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0,-wnum*two*al(5)*anu(5)) 

fftem=exp(ffteml ) 

ff(5)=(one+qq(5)*fftem)/(-one+qq(5)*fftem) 
anu(4)=sqrt(eps(4)-(sin(phi))*  *2) 
qqtem=(bh(5)/bh(4))*(anu(4)/anu(5))*ff(5) 
qq(4)=(qqtem+one)/(qqtem-one) 
ffteml=cmplx(0.0,-wnum*two*al(4)*anu(4)) 
fftem=exp(ffteml ) 

ff(4)=(one+qq(4)*fftem)/(-one+qq(4)*fftem) 

anu(3)=sqrt(eps(3)-(sm(phi))**2) 
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qqtein=(bh(4)/bh(3))*(anu(3)/anu(4))*ff(4) 
qq(3)=(qqtem+one)/(qqtem-one) 
ffteml=cmplx(0.0,-wnum*two*al(3)*anu(3)) 
fftem=exp(ffteml ) 

ff(3)=(one+qq(3)*fftem)/(-one+qq(3)*fftem) 

anu(2)=sqrt(eps(2)-(sm(phi))**2) 

qqtem=(bh(3)/bh(2))*(anu(2)/anu(3))*ff(3) 

qq(2)=(qqtem+one)/(qqtem-one) 

ffteinl=cmplx(0.0,-wnum*two*al(2)*anu(2)) 

fftem=exp(ffteml ) 

ff(2)=(one+qq(2)*fftem)/(-one+qq(2)*fftem) 

qqtem=(cos(phi)/anu(2))*ff(2) 

qq(l)=(qqtem+bh(l)/bh(2))/(qqtem-bh(l)/bh(2)) 

ffteml=cmplx(0.0,-wnuin*two*al(l)*cos(phi)) 
fftem=exp(ffteinl ) 

ff(l)=(one+qq(l)*fftem)/(-one+qq(l)*fftem) 

ref=(ff(l)+bhl/bh(l))/(ff(l)-bhl/bh(l)) 

refr=real(ref) 

refi=aimag(ref) 

refabs=sqrt(refr**2+refi**2) 

if  (refabs.GT.valuel)  then 
value  l=refabs 
end  if 
10  continue 
20  continue 
30  continue 
f=valuel 
return 
end 


Appendix  D 

COMPUTER  PROGRAM  “HE^C22” 


PROGRAM  DESCRIPTION 

“HINC22”  uses  the  optimum  impedance  matching  network  parameters  generated  by  “HOPT22”  for 
the  matching  network  of  Fig.  8(b),  which  uses  four  transformers.  “fflNC22”  uses  these  parameters  to 
compute  the  reflection  coefficient  at  each  frequency,  incidence  angle,  and  dielectric  constant  of  interest. 
The  program  writes  outputs  of  its  calculation  into  the  data  file  "aaaaa.dat".  The  output  data  file  consists  of 
four  columns:  the  first  colunon  lists  the  frequency  in  GHz,  the  second  column  gives  the  absolute  value  of 
the  reflection  coefficient,  the  third  column  is  the  absolute  value  of  the  reflection  coefficient  in  dB,  and  the 
fourth  column  is  the  phase  of  the  reflection  coefficient  in  radians.  Input  data  to  the  program  are 

a.  heights  of  the  waveguides  A,  B,  and  C; 

b.  length  of  the  waveguide  B  (length  of  the  ferroelectric);  waveguides  A  and  C  are  assumed  to  be 
semi-infinite; 

c.  dielectric  constant  of  the  ferroelectric; 

d.  incidence  (scan)  angle; 

e.  upper  and  lower  frequency  and  a  step  increment; 

f.  dielectric  constant  of  the  material  in  the  third  and  fourth  transformer; 

g.  lengths  of  all  four  transformers;  and 

h.  height  of  the  first  two  transformers. 

All  lengths  are  to  be  entered  in  inches,  frequencies  in  GHz,  and  angles  in  degrees. 

LISTING  OF  PROGRAM 
c  Symmetrical  transition:  4  transformers 


implicit  double  precision  (a-h,o-z) 

complex  qq,ff,ffteml,fftem,qqtem,ref 

dimension  qq(  1 00),ff(  1 00)  ,bh(  1 00)  ,al  ( 1 00) 

dimension  eps(  1 00),anu(  1 00),gam(  1 00) 

open  (unit=  1  ,file='aaaaa.dat',status— new') 

print  *,  "Results  are  in  file:  aaaaa.dat" 

one=1.0d0 

two=2.0d0 

ten=10.0d0 

api=acos(-one) 

api2=api/two 

apit=api 

api2t=api2 

constl=0.16933333d0 
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print  *,  "Enter  incidence  angle  in  degrees" 
read  (5,*)  phideg 
phi=phideg*api/l  SO.OdO 

print  *,  "Enter  low  frequency,  upper  frequency,  freq.step" 
read  (5,*)  freql,frequ,frstep 

print  *,  "Enter  larger  waveguide  height  (in  inches)" 
read  (5,*)  bhl 

print  *,  "Enter  smaller  waveguide  height  (in  inches)" 
read  (5,*)  bhs 

print  *,  "Enter  smaller  waveguide  length  (in  inches)" 

read(5,=^)allls  .  , 

print  *,"Enter  smaller  waveguide  dielectnc  constant  (relative) 

read  (5,*)  epss 
nst=4 

nstl=2*nst+l 

nst2=2*nst+2 

bh(nst+l)=bhs 

eps(nst+l)=epss 

al(nst+l)=allls 

print  *,  "Enter  height  of  first  transformer  (in  inches)" 
read  (5,*)  bh(l) 

print  *,  "Enter  height  of  second  transformer  (in  inches)" 

read  (5,*)  bh(2) 

bh(3)=bhs 

bh(4)=bhs 

bh(6)=bh(4) 

bh(7)=bh(3) 

bh(8)=bh(2) 

bh(9)=bh(l) 

print  *,  "Enter  length  of  first  transformer  (in  inches)” 
read  (5,*)  al(l) 

print  *,  "Enter  length  of  second  transformer  (in  inches)" 
read  (5,*)  al(2) 

print  *,  "Enter  length  of  third  transformer  (in  inches)" 
read  (5,*)  al(3) 

print  *,  "Enter  length  of  fourth  transformer  (in  inches)" 

read  (5,*)  al(4) 

al(6)=al(4) 

al(7)=al(3) 

al(8)=al(2) 

al(9)=al(l) 

print  *,  "Enter  dielectric  constant  of  third  transformer" 
read  (5.*)  eps(3) 

print  *,  "Enter  dielectric  constant  of  fourth  transformer" 

read  (5,*)  eps(4) 

eps(l)=1.0d0 

eps(2)=1.0d0 

eps(6)=eps(4) 

eps(7)=eps(3) 

eps(8)=eps(2) 
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eps(9)=eps(l) 

numfi^(frequ-freql)/frstep 

print  *,  'nunifi:^',nunifr 

do  10  ifn=0,numfr 

fr=freql+frstep*ifr 

wnum=api*constl  *fr 

cccc  1  =(bhl/bh(9)-one)/(bhl/bh(9)+one) 

qq(9)=cmplx(cccc  1 ,0.0d0) 

ffteml=cmpbc(0.0d0,-wnum*two*al(9)*dcos(phi)) 

fftem=exp(ffteml) 

ff(9)=(one+qq(9)*fftem)/(-one+qq(9)*ffteni) 
anu(8)=sqrt(eps(8)-(sin(phi))  *  *2) 
qqtem=(bh(9)/bh(8))*(anu(8)/cos(phi))*ff(9) 
qq(8)=(qqtem+one)/(qqtem-one) 
ffteml=cmplx(0.0d0,-wnum*two*al(8)*anu(8)) 
fftem=exp(ffteml ) 

ff(8)=(one+qq(8)*fftem)/(-one+qq(8)*ffiem) 

anu(7)=sqrt(eps(7)-(sin(phi))  *  *2) 

qqtem=(bh(8)/bh(7))*(anu(7)/anu(8))*ff(8) 

qq(7)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0d0,-wnum*two*al(7)*anu(7)) 

fftem=exp(ffteml ) 

ff(7)=(one+qq(7)  *fftem)/(-one+qq(7) 

anu(6)=sqrt(eps(6)-(sin(phi))**2) 

qqtem=(bh(7)/bh(6))*(anu(6)/anu(7))*ff(7) 

qq(6)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0d0,-wnum*two*al(6)*anu(6)) 

fftem=exp(ffteml ) 

ff(6)=(one+qq(6)  *fftem)/(-one+qq(6)  *fftem) 

anu(5)=sqrt(eps(5)-(sin(phi))**2) 

qqtem=(bh(6)/bh(5))*(anu(5)/anu(6))*ff(6) 

qq(5)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0d0,-wnum*two*al(5)*anu(5)) 

fftem=exp(ffteml) 

ff(5)=(one+qq(5)*f¥tem)/(-one+qq(5)*fftem) 

anu(4)=sqrt(eps(4)-(sin(phi))**2) 

qqtem=(bh(5)/bh(4))*(anu(4)/anu(5))*ff(5) 

qq(4)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(0.0d0,-wnum*two*al(4)*anu(4)) 

fftem=exp(ffteml ) 

ff(4)=(one+qq(4)  *fftem)/  (-one-i-qq(4)*f¥tem) 

anu(3)=sqrt(eps(3)-(sin(phi))**2) 

qqtem=(bh(4)/bh(3))*(anu(3)/anu(4))*ff(4) 

qq(3)=(qqtem+one)/(qqtem-one) 

ffteml=cmplx(O.OdO,-wnum*two*al(3)*anu(3)) 

fftem=exp(ffteml) 

ff(3)=(one+qq(3)*fftem)/(-one+qq(3)*fftem) 

anu(2)=sqrt(eps(2)-(sin(phi))**2) 

qqtem=(bh(3)/bh(2))*(anu(2)/anu(3))*ff(3) 

qq(2)=(qqtem+one)/(qqtem-one) 
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ffteml=cmplx(0.0d0,-wnum*two*al(2)*anu(2)) 
fftem=exp(ffteinl ) 

ff(2)=(one+qq(2)*ffteni)/(-one+qq(2)*fftem) 

qqtem=(cos(phi)/anu(2))*ff(2) 

qq(l)=(qqtem+bh(l)/bh(2))/(qqtem-bh(l)/bh(2)) 

ffteml=cmplx(O.OdO,-wnum*two*al(l)*cos(phi)) 

fftein=exp(ffteml ) 

ff(l)=(one+qq(l)*fftem)/(-one+qq(l)*fftem) 

ref=(ff(l)+bhl/bh(l))/(ff(l)-bhl/bh(l)) 

refr=real(ref) 

refi=aimag(ref) 

refabs=sqit(refr**2+refi**2) 

refdb=20.{)d0*log  1 0(refabs) 
refarg=atan(refi/refr) 
write  (1,20)  fr,refabs,refdb,refarg 
20  fomiat(lx,f5.2,lx,f8.5,lx,f6.2, 

$lx,f8.3) 

10  continue 
stop 
end 


